arXiv:1501.06532v2 [cond-mat.supr-con] 13 Mar 2015 


Electronic heat current rectification in hybrid superconducting devices 

Antonio Fornieri, 1, a) Marfa Jose Martfnez-Perez , 2 and Francesco Giazotto 1,b) 

l) NEST, Istituto Nanoscienze-CNR and Scuola Normale Superiore, 1-56127 Pisa, Italy 
21 Physikalisches Institut - Experimentalphysik II Universitat Tubingen, D-72076 Tubingen, 

Germany 

In this work, we review and expand recent theoretical proposals for the realization of electronic thermal diodes based 
on tunnel-junctions of normal metal and superconducting thin films. Starting from the basic rectifying properties of 
a single hybrid tunnel junction, we will show how the rectification efficiency can be largely increased by combining 
multiple junctions in an asymmetric chain of tunnel-coupled islands. We propose three different designs, analyzing 
their performance and their potential advantages. Besides being relevant from a fundamental physics point of view, this 
kind of devices might find important technological application as fundamental building blocks in solid-state thermal 
nanocircuits and in general-purpose cryogenic electronic applications requiring energy management. 


I. INTRODUCTION 

A thermal rectifier 1-3 can be defined as a device that con¬ 
nects asymmetrically two thermal reservoirs: the heat current 
transmitted through the diode depends on the sign of the tem¬ 
perature bias imposed to the reservoirs. This non-linear device 
represents the thermal counterpart of the well-known electric 
diode, which helped the extraordinary evolution of modern 
electronics together with the transistor. The realization of ef¬ 
ficient thermal rectifiers would represent a giant leap for the 
control of heat currents at the nanoscale, 4-6 boosting emerging 
fields such as coherent caloritronics, 7-9 nanophononics, 6 and 
thermal logic. 6 Furthermore, a great number of nanoscience 
fields, including solid state cooling, 4 ultrasensitive cryogenic 
radiation detection 4, 0 and quantum information, 11,12 might 
strongly benefit from the possibility of releasing the dissipated 
power to the thermal bath in an unidirectional way. 

A highly efficient thermal diode should provide differences 
of at least one order of magnitude between the heat cur¬ 
rent transmitted in the forward temperature-bias configura¬ 
tion, /f w , and that generated upon temperature bias reversal, 
/ rev . This is equivalent to say that the rectification efficiency 

( 1 ) 

Jrev 


must be 1 or <C 1. 

Since the first theoretical proposal, 2 several groups have 
put a great effort into envisioning different designs for ther¬ 
mal rectifiers dealing with phonons, 13-16 electrons 17-22 and 
photons. 23 Alongside these theoretical works, promising ex¬ 
perimental results were obtained in systems that exploited 
phononic 24-26 and electronic 27 thermal currents. Neverthe¬ 
less, a maximum 51 ^ 1.4 has been reported 25 and more effi¬ 
cient rectification mechanisms are required to realize compet¬ 
itive thermal diodes. 

In this article, we will focus on the rectification perfor¬ 
mance of devices based on tunnel-junctions between normal 
metal (N) and superconducting (S) thin films at low temper¬ 
atures. First, we will review the intrinsic properties of single 


^Electronic mail: antonio.fornieri@sns.it 
b) Electronic mail: francesco.giazotto@sns.it 


NIN, NIS and S 1 IS 2 junctions 28,29 (where I stands for an in¬ 
sulating barrier and Si and S 2 represent two different super¬ 
conducting electrodes). Then, we will analyze possible im¬ 
provements of the rectification efficiency when these struc¬ 
tures are combined together forming an asymmetric chain of 
tunnel-coupled islands. 30,31 In particular, as we shall argue, an 
asymmetric coupling to the phonon bath can strongly enhance 
the efficiency of the device, providing outstanding values of 
51 with realistic parameters. Very recently, the latter mecha¬ 
nism has been experimentally demonstrated in a hybrid device 
exhibiting a maximum rectification efficiency of ~ 140. 31 


II. NIN AND NIS JUNCTIONS 


We shall start, first of all, by describing the equations 
governing heat transport in simple NIN and NIS junctions. 
In these opening sections, we focus on the electronic heat 
currents only, and neglect possible contributions by lattice 
phonons. The latter will be discussed in the analysis of chain 
diodes. 

If we consider two N electrodes residing at electronic tem¬ 
peratures 7h 0 t and T co \& (with 7h 0 t > T co \ d) coupled by means 
of a tunnel junction, the stationary electronic thermal current 
flowing through the junction can be written as: 28 


^NIN ( ^hot 5 ^cold ) 


6e 2 R T [hoi 



( 2 ) 


where Ry is the contact resistance, e is the electron charge and 
k b is the Boltzmann’s constant. Equation (2) clearly shows 
that no rectification is possible in a full normal-metal tunnel 
junction, since |/NIN(^hot 5 ^cold)| = I^NIN^cold, 7hot)|- 

However, if we substitute the second N electrode with a 
superconductor, the thermal current flowing through the NIS 
junction becomes: 4,32 


^NIS (^hot? ^cold) — Jq dEEyi(E, 7c 0 id) 


( 3 ) 


x\f(E,T hot )-f(E,T cold )\. 

Here, N(£,r) = 91 \e + iT /y/(E + iT) 2 - A 2 (T)j is 

the smeared (by non-zero T) normalized Bardeen- 
Cooper-Schrieffer (BCS) density of states (DOS) in the 
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FIG. 1. (a) Thermal current J vs. T^ ot flowing through a NIN junction (dashed line) and through a NIS junction (solid lines). The curves are 
calculated for Rj = 1 k£2, T co ^ = 0.01 T c and T c = 1.4 K. The inset highlights different trends of J in the same cases as those displayed in the 
main panel for small values of T^ ot on a logarithmic scale, (b) Rectification efficiency 31 vs. T^ot of a NIS junction for different values of T co 
Inset: optimal rectification efficiency $ 0 pt vs. T co \ ( Horizontal dashed lines indicate 31 = 1. 


superconductor, 33 35 A (T) is the temperature-dependent 
superconducting energy gap with a critical temperature 

T c = A(0)/(1.764fe B ) and f(E,T) = [1 + exp(^)]- 1 is 

the Fermi-Dirac distribution function. In the following, 
we will assume T ~ 10 _4 A(0), which describes realistic 
NIS junctions. 31,34,35 We also define /f w = |/Nis(^hot, 7ioid)| 
and / rev = |*/nis (^coicb ^hot) I for the forward and the reverse 
configuration, respectively. Figure 1(a) shows the influence 
of the superconducting DOS on J: at a first glance, /f w imme¬ 
diately appears different from / rev , thanks to the temperature 
dependence of A (T), which breaks the thermal symmetry 
of the tunnel junction. For 7h ot > 0.4!T C the value of A (T) 
in the reverse configuration starts to decrease and i rev gets 
significantly larger than Jf w , leading to 31 < 1, as shown in 
Fig. 1(b). On the contrary, if 7h 0t is raised above T c , heat flux 
from N to S becomes preferred, resulting in 31 > 1. As T co \& 
is increased, the reverse rectification regime gets gradually 
suppressed, while the forward regime becomes more efficient. 
We can highlight this behavior by introducing the optimal 
rectification efficiency lR opt , defined as that corresponding to 
the maximum value between 31 and 1/31 at a given r co id. This 
quantity provides the best working point of the diode, which 
can occur in the forward (if $ opt > 1) or in the reverse (if 
^opt < 1) regime of rectification. As displayed in the inset 
of Fig. 1(b), for Tcoid < 0.5T C the reverse regime is the most 
efficient, while at higher temperatures the forward regime is 
favored. Finally, at T co \^ = T c the system turns into a NIN 
junction and no rectification is possible. 

It is also worth noting that at low values of 7h ot and r co id, 
/nis is suppressed by a factor T with respect to /nin and does 


not depend on the sign of the temperature bias [see inset of 
Fig. 1(a)], indicating the effectiveness of A (T) as a thermal 
bottleneck. On the other hand, when 7hot and T co ^ increase, 
/(£,7hot) — f{E,T C0 \d) becomes significantly different from 
zero at energies higher than the smeared A (T): this corre¬ 
sponds to a rapid increase of both /f w and / rev , as indicated by 
the arrow in the inset. Thus, the thermal bottleneck effect gen¬ 
erated by the superconducting energy gap decays quickly by 
raising temperatures. Although this change of behavior is not 
relevant for the rectifying properties of a single NIS junction, 
it will be important to understand the operation of a NINISIN 
chain. 


III. S!lS 2 JUNCTIONS 

We now move to the case in which both the tunnel-coupled 
electrodes are superconductors, i.e., a SiIS 2 junction. 29 The 
latter can be thermally biased by setting the quasiparticle tem¬ 
perature of one electrode to 7h 0 t and that of the other to T co \^. 
Moreover, we assume that Si and S 2 are characterized by 
different superconducting energy gaps Ai and A 2 and by a 
phase difference (p = (p 2 — <Pi. Then, the forward and re¬ 
verse total heat currents flowing through the Josephson junc¬ 
tion read: 32,36 

^SiIS 2 ,fw(rev) = Jqp [^ot(cold)? ^cold(hot)] 

— Ant[^hot(cold)5 ^cold(hot)] cos <P- 

Here, the term / qp accounts for the energy carried by 
quasiparticles 32 ’ 37-41 and is equivalent to Eq. 3 for the present 
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FIG. 2. (a) Heat currents /f w and / rev vs. T^ot flowing through a S 1 IS 2 junction for <p = 0, n. The curves are calculated for 8 = 0.75, Rj = 1 
k£2, T c oid = 0.01 T c \ and T c \ = 1.4 K. (b) Contour plot showing R as a function of 8 and T^ot for (jp = n and / co id = 0.01 T c \. (c) Three selected 
profiles of R vs. 8 for the values of T^ot corresponding to the dashed lines in panel (b). (d) R vs. T^ot for 8 = 0.75 and for different values of 
<P at ^cold = 0.01 T c \. The horizontal dashed line indicates R= 1. 


system: 


^qp(^hoO ^cold)— 2d [ dEEJ^\ [E 2 {R ffcold) _ 

e z Rj Jo (5) 

x[/(7h 0t )-/(r C0ld )], 

where 3Sfi and N 2 are the smeared normalized superconduct¬ 
ing DOS of Si and S 2 , respectively. On the other hand, 7 mt is 
the phase-dependent component of the heat current: 32,38 ^ 1 

j int ( ^hot 5 ^cold ) = 2 d / dEEJA \ (E, T\ lo \)JA2(E ^ T co \&) 

e z Rj Jo (6) 

x \f(j hot) -/(^cold)]. 


where M k {EJ) = 3 -iA k (T) / J (E+ iT k ) 2 - A 2 (T) 


is 


the Cooper pair BCS density of states in S^ at temperature 
71, 42 ~ 10 _4 A^(0) and k = 1,2. This term is peculiar to 
the Josephson effect (so it vanishes if one or both of the su¬ 
perconducting energy gaps are null) and originates from tun¬ 
neling processes through the junctions involving both Cooper 
pairs and quasiparticles. 32,38 Depending on <p, J mi can flow in 
opposite direction with respect to that imposed by the thermal 
gradient, but the total heat current still flows from the hot to 
the cold reservoir, preserving the second principle of thermo¬ 
dynamics. This was experimentally demonstrated in Ref. 7. 

In the following, we shall assume, for clarity, that 8 = 
A 2 /A 1 < 1. Figure 2(a) shows the behavior of Jf w and J rev 
vs. Thot for r co i d = 0.01 T c i (T c 1 being the critical tempera¬ 
ture of Si) and 8 = 0.75. The difference between the heat 
currents is particularly evident for (p = tc\ in the forward con¬ 
figuration, /f w exhibits a sharp peak at Thot — 0.777^1, due to 


the matching of singularities in superconducting DOSs when 
Ai(7hot) = A 2 (r cold ). At higher values of 7h ot , Ai(7h ot ) < 
A 2 (r co i d ) and the energy transmission through the junction 
is reduced, leading to a negative thermal differential conduc¬ 
tance region, which is further enhanced by the gradual sup¬ 
pression of /int as Thot approaches T c \. In the reverse con¬ 
figuration, instead, / rev presents just one cusp caused by the 
vanishing of J mi when Thot = T C 2 = 0.75/ c i. Then, the rec¬ 
tification efficiency is maximum when the aforementioned 
points are perfectly aligned, i.e., R 0 pt — 7.5 for 8 = 0.75 and 
Th 0t — 0.77/d at r co i d = 0.01 T c i, as shown in Figs. 2(b) and 
2(c). As r co id is increased, R opt decreases monotonically to 
~ 1.3 at r co id = 0.7/ c i, then it becomes < 1 and finally reaches 
the identity for T co id = T c \ (data not shown). More interest¬ 
ingly, the rectification efficiency can be finely tuned by con¬ 
trolling the phase bias across the junction. As a matter of fact, 
if (p = 0, R is not only reduced by one order of magnitude but 
it also inverts the favored rectification regime, as displayed by 
Fig. 2(d). This effect can be readily understood by observing 
the trend of Jf w and J rev vs. 7h ot for (p = 0 [see Fig. 2(a)]. In 
this case, the negative sign in front of J[ nt in Eq. 6 perfectly 
cancels the singularity-matching peak in the forward configu¬ 
ration. 

To conclude the analysis of the intrinsic properties of NIS 
and S 1 IS 2 junctions, it is worthwhile to summarize the con¬ 
ditions necessary to achieve heat rectification: (1) two differ¬ 
ent DOSs must be tunnel coupled (i.e., Ai / A 2 ) and (2) at 
least one of them must be strongly temperature dependent in 
the range of operation. Furthermore, a fairly large tempera¬ 
ture gradient is needed, since heat rectification is absent in the 
linear response regime, i.e., when 8T = 7h ot — ^coid = 










































4 


(T bol + T mli )/ 2 29 


IV. NININ CHAIN 

The rectifying performance of the discussed tunnel junc¬ 
tions can be largely improved if they are combined in an asym¬ 
metric chain of tunnel-coupled electrodes that can exchange 
energy with an independent phonon bath residing at tempera¬ 
ture 7b a th. We stress that we are concerned with the heat cur¬ 
rent carried by electrons only. We assume that lattice phonons 
present in the entire structure are thermalized, with substrate 
phonons residing at r bat h, and are therefore not responsible for 
any heat transport. This assumption is expected to hold be¬ 
cause the Kapitza resistance between thin metallic films and 
the substrate is vanishingly small at low temperatures . 7,8,31,43 

We shall begin with a simple N 1 IN 2 IN 3 chain , 30 which, as 
we will show, represents a fundamental element to understand 
the working principles of chain diodes. Ni and N 3 act as ther¬ 
mal reservoirs and are used to impose a temperature gradient 
across the device. Since is defined under the condition of 
equal temperature bias, these electrodes must be identical and 
equally coupled to the phonon bath. Thus, 7b at h now plays 
the role embodied by T co \^ in the single junction analysis. N 2 
is connected to Ni and N 3 by means of two tunnel junctions 
characterized by resistances R\ and R 2 , respectively. For sim¬ 
plicity, we assume R\ equal to 1 k Q and consider only the 
parameter r = R 2 /R 1 > 1, accounting for the asymmetry of 
the chain. N 2 is the core of the diode and controls the heat 
flow from one reservoir to the other by releasing energy to 
the phonon bath. In the forward configuration, the electronic 
temperature of Ni is set to 7h 0t > 7b a th, leading to electronic 
thermal currents (see Eq. 2) / 2 ,fw and /f w flowing into N 2 and 
N 3 , respectively. On the contrary, in the reverse configuration 
the electronic temperature of N 3 is set to 7h 0 t, generating heat 
currents / 2 ,rev and / rev flowing into N 2 and Ni, respectively. 

Besides electronic thermal currents, we must take into ac¬ 
count the heat exchanged by electrons in the metal with lattice 
phonons : 8,44 


/ N ,ph(7\ Tbath) = EV(r n - r b " ath ). (7) 

Here £ is the material-dependent quasiparticle-phonon cou¬ 
pling constant, V is the volume of the electrode and n is 
the characteristic exponent of the material. In this work we 
will consider two materials that are commonly exploited to 
realize N electrodes in nanostructures, i.e., copper (Cu) and 
manganese-doped aluminum (AIMn). The former is typically 
characterized by Ecu = 3 x 10 9 WK - 5 m -3 and nc u = 5 ,' 4,7 
while the latter exhibits £aimii = 4 x 10 9 WK - 6 m -3 and 
^ AiMn = 6 . 8,3 1,44 Furthermore, we assume that all the elec¬ 
trodes of the chain have a volume Vn = 1 x 10“ 20 m- 3 . 
Equations 2 and 7 can be used to elaborate a thermal model 
accounting for heat transport through the device. The model 
is sketched in Fig. 3(a) and describes the forward tem¬ 
perature bias configuration, in which the electrodes of the 
chain reside at temperatures 7hot > 72 ,f w > 7f w > 7t> a th- Here, 
72 , fw and 7} w represent the electronic temperatures of N 2 
and N 3 , respectively. The terms 72,fw = ^NiN(7hot, ^ 2 ,fw) and 


Jfw = Tnin (^2,fw ? ^fw) account for the heat transferred from 
Ni to N 3 . The reservoirs can release energy to the phonon 
bath by means of /i, ou t and / 3 , 0 m• Photon-mediated thermal 
transport , 45-47 owing to poor impedence matching, as well as 
pure phononic heat currents are neglected in our analysis . 7,8,31 
We can now write a system of energy-balance equations that 
account for the detailed thermal budget in N 2 and N 3 by set¬ 
ting to zero the sum of all the incoming and outgoing heat 
currents: 

^2,fw(^hot5 ^2,fw) — Aw(^2,fw? ^fw) — *^ 2 ,out(^ 2 ,fw 5 ^bath) = 0 

^fw(^2,fw^fw) —/3,out(7fwj^bath) = 0- (8) 

Here / 2 ,out(^ 2 ,fw 7 ^bath) is the heat current that flows from 
N 2 to the phonon bath. In Eqs. 8 we can set 7h ot and r batb 
as independent variables and calculate the resulting ? 2 ,f w and 
7f w . Another system of energy-balance equations can be writ¬ 
ten and solved for the reverse configuration , 48 in which N 2 and 
Ni reach electronic temperatures 72,rev and T rev , respectively. 
Finally, we can extract the values of 7f w and 7 rev , thereby ob¬ 
taining R. 

In order to grasp the essential physics underlying this de¬ 
vice, it is instructive to consider an ideal case, in which Ni, 
N 2 and N 3 can release energy to the phonon bath only through 
three tunnel-coupled N electrodes, acting as thermalizing cold 
fingers 30 and labeled as Fi, F 2 and F 3 , respectively. We as¬ 
sume that electrons in Fi, F 2 and F 3 are strongly coupled to 
the phonon bath, so they reside at T bat h. As we shall argue, 
F 2 is fundamental to achieve a high rectification performance 
and is characterized by a resistance Rp, 2 , while Fi and F 3 
must have the same tunnel junction resistance Rp,i- We ig¬ 
nore at this stage contributions due to the direct quasiparticle- 
phonon coupling (see Eq. 7) in Ni, N 2 and N 3 , in or¬ 
der to obtain simple analytic results for the rectification effi¬ 
ciency. As a matter of fact, this ideal device can be easily de¬ 
scribed by substituting ^out^.fw, ^bath) = -^NIN^fw, ^bath) 
and 1 / 3 ,out(?fw 7 Tbath) = *^NiN(^fw: ^bath) hi Eqs. 8 . Remarkably, 
the resulting expression for the rectification efficiency does 
not depend on 7h 0t nor T bat h, but only on the resistances of the 
tunnel junctions defining our NININ chain: 

gj _ Rl (^2 +^F, 2 ) +^ 2 (^F ,1 +^F, 2 ) +^F, 1 ^F ,2 
^F, 2(^2 +7?F,l) +^l (^2 +*F,1 +^F, 2 ) 

Equation 9 is particularly transparent in two limit cases: first, 
if Rp ,2 —>■ 00 , ft 1 for every value of r. This clearly demon¬ 
strates that no rectification is achievable if the central island 
of the chain is not coupled to the phonon bath, regardless the 
asymmetry of the device. This result holds true for a general 
design of a NININ chain, provided that Ni and N 3 are identi¬ 
cally coupled to the phonon bath . 30 

On the other hand, if Rp,i —^ 00 , Eq. 9 is simplified as fol¬ 
lows: 


rt> _^2+^F,2 

^Mio-ph — „ . „ • 

K\ +/vp ,2 


( 10 ) 


This expression summarizes the conditions needed to obtain 
an highly-efficient NININ diode: (1) the thermal symmetry of 
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FIG. 3. (a) Thermal model for the NININ chain in the forward thermal bias configuration, which is characterized by T^ot > ^2,fw > 7f w > T batb . 
The most efficient design for the diode consists in tunnel-coupling a cold finger F 2 to N 2 and in minimizing the effect of the quasiparticle- 
phonon coupling in all the electrodes composing the chain (see text), (b) 01 vs. 7h ot for different values of the cold finger resistance Rp ,2 at 
T^ath = 50 mK. (c) Diode’s output temperatures 7} w (solid lines) and r rev (dashed lines) vs. T^ot at different values of r batb . (d) 0l opt vs. r batb 
for two different N materials. In panel (c) and (d) we set 2 = 500 £1 All the results have been obtained for r = 100. 


the device must be broken (r 1 ) and ( 2 ) N 2 must be very 
well coupled to the phonon bath with respect to Ni and N 3 . 
A natural trade off arising from the fulfilment of these condi¬ 
tions is the reduction of the thermal efficiency of the device 
I? = */f w /^ 2 ,fw < 1 , he., the fraction of energy that is actually 
transferred from one reservoir to the other in the transmissive 
regime . 30 Nevertheless, the device parameters can be designed 
in order to maximize the required performance in terms of the 
global efficiency Otrj , 30 which can range from 0 to 00 , since 
01 has no upper limit. Finally, this ideal case demonstrates 
that such a thermal rectifier can operate efficiently also in the 
regime for T hot -A r bath . 

In a more realistic case, every electrode in the chain can 
exchange energy with the phonon bath by means of the 
quasiparticle-phonon coupling, which introduces a tempera¬ 
ture dependence in the rectification efficiency. It is then useful 
to rewrite 01 1 as a general condition for temperatures: 


3= _ ,wS / ,w >1, 

rT rev 8T rev 


(ID 


where <5r fw ( rev ) = r 2?fw ( rev ) - T fw ( rev ) and the mean tempera¬ 
tures T fw ( rev) = (r 2?fw ( rev ) + 7fw(rev))/2. In order to maximize 
01, in the reverse configuration the electronic temperatures of 
Ni and N 2 must be similar and close to the lowest temperature 
in the system, i.e., r bat h- This can be easily done by setting 
r > 1 and by coupling N 2 to the phonon bath. On the other 


hand, the coupling between N 2 and the phonon bath should 
have a limited impact on heat transport in the forward con¬ 
figuration. From Eq. 7, it appears evident that /n ?p 1 i is n °t 
well-suited to satisfy the aforementioned requirements, since 
it is not much effective at low temperatures, while it becomes 
strongly invasive at high temperatures. Thus, it turns out that 
the most efficient design for a NININ chain diode is obtained 
by coupling just one thermalizing cold finger F 2 to N 2 and by 
minimizing the impact of the quasiparticle-phonon coupling 
in all the electrodes of the chain . 30 

Then, we can set /2,out(^2,fw5 ^bath) = ^N,ph(^2,fw5 ^bath) H - 
*^NIN(^2,fw5 -^bath) and /3 ? out(^fw5 ^bath) = TN ?p h(7f w , -^bath) hi the 

energy-balance equations. Figure 3(b) shows 01 as a func¬ 
tion of Thot for different values of Rp ? 2 at T bat h = 50 mK. Re¬ 
markably, an optimal rectification efficiency 49 0i opt ~ 3000 
is obtained for Rp ,2 = 500 Q and r = 100. This configu¬ 
ration generates a maximum global efficiency pR7j] ma x — 8 , 
which can be optimized up to a value of about 9.6 by setting 
Rf ,2 = 1.2 k£2. We can also notice that for Thot -A r bat h, 01 
approaches IR no _ p h, which represents an asymptotic value for 
the rectification efficiency at low temperatures, i.e., when the 
quasiparticle-phonon coupling is almost uneffective. 

Figure 3(c) displays the rectifier’s output temperatures 7} w 
and T wv vs. 7h 0t for r = 100 and R? = 500 Q at different val¬ 
ues of r bath . The results point out a maximum difference of 
~ 200 mK between the forward and reverse configurations at 
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Tbath = 50 mK. Additionally, the behavior of output tempera¬ 
tures clearly highlights the effectiveness of F 2 in keeping r rev 
close to 7b at h, while /N,ph strongly limits the increasing rate of 
7f w as Thot and Tb a th become larger. 

Finally, Fig. 3(d) confirms the negative effect of the 
quasiparticle-phonon coupling on the performance of the 
diode. As a matter of fact, the behavior of R opt as a func¬ 
tion of Tbath shows that the T 5 dependence of /N,ph in Cu is 
extremely detrimental in the forward configuration and can 
reduce the rectification efficiency up to a factor of 10 . 


V. NINISIN CHAIN 


Despite its conceptual simplicity and potential insensitivity 
to magnetic fields, the NININ chain diode might be demand¬ 
ing from a fabrication point of view because of the required 
high asymmetry between R\ and R 2 . A good alternative is 
represented by a N 1 IN 2 ISIN 3 chain, in which all the tunnel- 
junction normal-state resistances R\, R 2 and R 3 are equal and 
the thermal symmetry is broken by the presence of the energy 
gap in the superconducting density of states. At low temper¬ 
atures, the latter plays a role equivalent to the one embodied 
by the largest tunnel-junction resistance in a NININ chain. In 
the following analysis, we will consider an S electrode made 
of aluminum with a critical temperature T c = 1.4 K. 

Emulating the most efficient design for a NININ chain, also 
in this case the N 2 electrode can be coupled to the phonon 
bath by means of a cold finger F residing at Tbath- On the 
other hand, it is worthwhile to note that a simple NISIN 
chain would not represent an efficient design because A (T) 
would prevent the central electrode to release energy to F 
in an efficient way, especially at low temperatures [see Fig. 
1(a)]. Figure 4(a) sketches the thermal model accounting for 
thermal transport through the device in the forward config¬ 
uration, characterized by T hot > r N?fw > Ts,fw > Tf w > r ba th- 
Here, r N _f w and T S f w label the electronic temperatures of 
N 2 and S, respectively. We assume that all the electrodes 
have the same volume V = 1 x 10 -20 m 3 and R\ = R 2 = 
R 3 = 1 kQ. Heat is transferred from Ni to N 3 by means 
of the heat currents / N iN(Thot,T^fw), /Nis(Tk,fw,Ts,f w ) and 
Jf w = *Tnis (Ts,fw? Tf w ). Electrode N 2 can release energy to the 
phonon bath via / 2 ,out = ^NiN(TN,fw? Tb a th) H - -TN,ph(TN ? fw? Tbath)? 
while N 3 is coupled to Tb at h through /N,ph(Tf w , Tbath)- Finally, 
S is affected by the quasiparticle-phonon coupling as well, 
which can be written as : 50 


-/s,ph(T\ Tbath) — 


LV 


x Le,e+£ 


coth 


96C(5)*4 J- 
£ 


L dEE L 


dee 2 sgn(e) 


2kg Tb a th 


f AL) _ \ _ 1 1 

\JE JE+e) JE JE+e I" 1 


( 12 ) 


where fjp = f(—E,T) — f(E,T) and the phonons are 
assumed to be in thermal equilibrium with occupation 

n(e, Tb a th) = [exp(j^r^) - 1] _1 . The factor L E . )E ' = 


N(E,r)3sr(E',r) 


1 - 


A 2 (T) 

EE' 


. In our case, we assume Lai = 


0.3 x 10 9 WK -5 m 3 . 4 With all these ingredients we can write 
and solve the energy balance equations for the forward and re¬ 
verse configuration, thus allowing us to analyze the operation 
of the diode vs. Th ot and Tbath for different parameters. 

Figure 4(b) displays the behavior of Jl as a function Th ot 
for different values of the resistance Re between F and N 2 
at Tbath = 50 mK. At low values of Th 0t the rectification ef¬ 
ficiency can reach values as high as 5000 for Re = 0.5 k£2, 
but for increasing temperatures it dramatically drops and ap¬ 
proaches unity for Th 0t > 1 K. This trend can be easily under¬ 
stood by first focusing on the case of Re > 10 kfZ The upper 
panel of Fig. 4(c) shows the diode’s output currents /f w ( rev ) 
vs. Th ot for Rf = 50 kfl and directly explains the Tl curve 
(solid line) reported in the lower panel. As already noticed 
at the end of Sect. II, for sufficiently high thermal bias an 
S electrode loses its effectiveness as a thermal bottleneck, as 
can be noticed from the change in the derivative of the out¬ 
put currents for Thot > 0.2 K in both the forward and reverse 
configurations. This change of behavior produces a maximum 
in the rectification efficiency at Th 0t ~ 0.15 K, after which R 
rapidly decreases. Then, the first part of the curve can be com¬ 
pared to the efficiency (dashed line) of an equivalent NINININ 
chain with identical volumes of electrodes and a resistance 
asymmetry r = R 3 /R 1 = 12000, which is about the value of 
A(0)/r = 10 4 . On the other hand, for Re < 10 k£2, Ts, rcv can 
be significantly higher than T$ f w at the same value of Th 0t . In 
particular, for Re = 0.5 kfl in the reverse configuration S loses 
its properties of thermal bottleneck for Th ot > 0.15 K, while in 
the forward configuration A blocks efficiently the thermal flux 
up to Th ot — 0-27 K. This discrepancy in the "breaking points" 
of the S electrode is reflected in the behavior of % which ex¬ 
hibits a shoulder after the first maximum [see Fig.4(b)]. Fur¬ 
thermore, it is worth noting that for Re > 2.5 k£2 the device 
can operate also in the reverse regime of rectification at suffi¬ 
ciently high Thot, reaching values of ‘R ^ 0.8, as expected for 
a simple NIS junction (see Sect. II). 

As Tbath is increased, the maximum value of R in the for¬ 
ward regime becomes strongly suppressed due to the invasive 
enhancement of the quasiparticle-phonon coupling. This can 
be indirectly observed in Fig. 4(d), which shows the diode’s 
output temperatures T fw ( rev ) vs. Th ot for Re = 0.5 kD at three 
selected values of Tbath- At Tbath = 50 mK we obtain a max¬ 
imum temperature difference ST = Tf w — T rev of about 110 
mK, while at higher Tb a th it rapidly decreases in conjunction 
with the derivative of T fw ( rev ). 

Remarkably, the global efficiency of this device at Tbath = 
50 mK is comparable to the NININ case and can be optimized 
to [fRT?] max — 7.5 by setting Re = 7.5 kfZ 

Finally, it is also possible to use a different superconductor 
with a larger A (such as vanadium or niobium) to implement 
the S electrode of the chain. This can significantly increase 
the range of Th ot where the superconductor can act as a large 
thermal resistance, thus providing % 1. Nevertheless, the 

performance of the device vs. Tbath will not be significantly 
improved, since it is mainly determined by the disruptive ef¬ 
fect of the quasiparticle-phonon coupling in N electrodes. 

The effectiveness of this kind of device has been recently 
proven in the experiment reported in Ref. 31 . 
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FIG. 4. (a) Thermal model for the NINISIN chain in the forward thermal bias configuration, characterized by T^ot > ^N,fw > ^S,fw > 7f w > ^bath- 
(b) vs. Thot f° r different values of the cold finger resistance R p at T^ath = 50 mK. (c) The upper panel shows Jf w (solid line) and / rev (dashed 
line) vs. T^ot for the NINISIN chain. The corresponding is represented by the solid line in the lower panel, while the dashed line stands for 
the rectification efficiency of a NINININ chain with a resistance asymmetry r = 12000. In both the panels we set R p = 50 k£2 and T b at h = 50 
mK. (d) Diode’s output temperatures 7f w (solid lines) and r rev (dashed lines) vs. T^ot at different values of T^ath for = 0.5 


VI. NINIS2IS1IN chain 

In previous works , 9,29 we analyzed theoretically the per¬ 
formance of a SiIS 2 junction tunnel-coupled to two N elec¬ 
trodes acting as thermal reservoirs. This thermal diode ex¬ 
hibits a remarkable phase-modulation of the rectification ef¬ 
ficiency for T c 1 = 1.4 K and <5 = 0.75 at 7b at h = 10 mK: the 
best performance is reached at 7h 0 t = 1-15, where Jl can be 
varied from a maximum of ^ 4.4 (if cp = n) to a minimum 
of ~ 0.7 (if (p = 0). 29 Thus, the behavior of a NISiIS 2 IN 
chain can be fully tuned by varying the phase difference be¬ 
tween Si and S2, allowing to invert the direction of the rec¬ 
tification regime. Nevertheless, this device has some exper¬ 
imental drawbacks, since its optimal working point requires 
a large bias temperature, close to the critical temperature of 
aluminum. As a consequence, superconducting thermometers 
used to detect 7h ot , 7} w and r rev (as in Refs. 7,8,31) should be 
realized with a superconductor S3, characterized by an energy 
gap A3 ^ Ai,A 2 , thus complicating the experimental design 
of the rectifier. Therefore, we can think of employing an addi¬ 
tional N electrode connected to the phonon bath by means of 
a cold finger in order to improve IR opt and to move the optimal 
working point of the device towards lower values of Thot- As 
we will argue, the proposed design offers the possibility to in¬ 


crease $ 0 pt by more than 3 orders of magnitude with respect to 
the NISiIS 2 IN chain diode and to demonstrate experimentally 
the phase modulation of Ji in a relatively simpler way. How¬ 
ever, it prevents us inverting the rectification regime through 
the phase modulation. 

We consider a NiIN 2 IS 2 ISiIN 3 chain, in which N 2 is 
tunnel-coupled to a cold finger F residing at 7b at h, Si has a 
critical temperature T c 1 = 1.4 K and <5 = A 2 /Ai = 0.75. This 
design results to be the one which maximizes the effect of a <j 0 
variation on the rectification efficiency. Figure 5(a) displays 
the thermal model accounting for thermal transport through 
the device in the forward configuration, characterized by 
^hot ^ ^N,fw ^ 7s 2 ,fw ^ T s j f w > 7f w > Here, 7jsr,f w , ^Si,fw 
and 7s 2 ? f w label the electronic temperatures of N 2 , Si and S 2 , 
respectively. We assume that all the electrodes have the same 
volume V = 1 x 10 _2 ° m -3 and all the tunnel junction have 
a normal-state resistance R\ = R 2 = R 2 = R 4 = l kfL Heat 
is transferred from Ni to N 3 by means of the heat currents 
^nin(^ hot 1 ^N,fw)? ^Nis(^N,fwTs 2 ,fw)> ^S]is 2 (7s 2 ,fw,^s,,fw) and 
Jfw = Tnis(^S i,fw 5 7f w ). Electrode N 2 can release energy to the 
phonon bath via / 2 out = *^NIN(^N,fw? ^bath) H“ JN,ph{JN ,fw? ^bath)> 
while N 3 is coupled to 7b at h through /N,ph(7fw, ^bath)- Finally, 
Si and S 2 are affected by the quasiparticle-phonon coupling 
/ Si, 2 ,ph( 7 s l 52 fwTbath). In general, the junction S 2 ISi can be 
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FIG. 5. (a) Thermal model for the NINIS 2 IS 1 IN chain in the forward thermal bias configuration, characterized by 7h ot > > ^s 2 ,fw > 

7si f w > 7f w > Tbath- (b) R vs. 7h 0 t for different values of the cold finger resistance R p at T^ath = 50 mK. Solid lines correspond to the case in 
which (p = 0, while dashed lines represent R when (p = n. (c) Diode’s output temperatures 7f w (solid lines) and r rev (dashed lines) vs. T^ot for 
R? = 2.5 kD and for (p = 0, n at at h = 50 mK. (d) 5T = 7f w — r rev vs. (p for different values of 7h ot at ^bath = 50 mK and for R p = 2.5 k£2. 


phase biased through supercurrent injection or by applying an 
external magnetic flux . 29 In order to achieve a complete mod¬ 
ulation of R, (p must vary between 0 and n. This can be ob¬ 
tained by using a radio frequency superconducting quantum 
interference device (rf SQUID), but this would again require 
the use of a clean contact between the S 2 IS 1 junction and a 
third superconductor S 3 (with A 3 Ai,A 2 ) in order to sup¬ 
press heat losses . 9,29 A good alternative is represented by an 
asymmetric direct current SQUID formed by different Joseph- 
son junctions (JJs), of which one corresponds to the S 2 IS 1 
junction. The latter can be biased up to (p = n if its Josephson 
inductance 51 is the largest in the SQUID, i.e., its Josephson 
critical current is much lower than those of the other JJs. 

As done in previous cases, we can now write and solve the 
energy-balance equations accounting for the thermal budget in 
every electrode of the chain. Thus, we can obtain the behavior 
of the thermal diode vs. 7h ot , 7b at h and (p for different values 
of the parameters. 

Figure 5(b) shows R vs. 7h 0 t for different values of the 
cold finger resistance R f at 7b at h = 50 mK. Solid lines repre¬ 
sent the rectification efficiency for (p = 0 , while dashed lines 
stand for that obtained in the case of (p = n. The values of R 
provided here are comparable to those obtained with the NIN- 
ISIN chain, but the phase dependence of /siis 2 represents an 
additional ingredient that strongly influences the performance 


of the diode over a large part of the 7h 0 t range. In particu¬ 
lar, for = 0.5 k Q we can reach a maximum relative vari¬ 
ation [$(0) — R(7r)\/R(0) ~ 77% at 7h ot = 0.13 K. From the 
point of view of output temperatures, though, the phase tun¬ 
ing of the device is clearly detectable only for 7h 0 t > 0.2 K. 
The effect of the phase biasing is especially visible on r rev , 
which is affected by a maximum variation of ~ 60 mK at 
Thot — 0.6 K. At this temperature also the phase modulation of 
8 T = 7J W — r rev exhibits the largest amplitude (of ~ 30 mK), 
as shown in Fig. 5(d). 

We finally notice that the global efficiency of a 
N 1 IN 2 IS 2 IS 1 IN 3 chain is comparable to those of other chain 
diodes, with a maximum pRrj] max ^ 6 for R? = 8 k £1. 


VII. SUMMARY AND FINAL REMARKS 

In summary, we first reviewed the intrinsic properties of su¬ 
perconducting hybrid junctions as thermal rectifiers. NIS and 
S 1 IS 2 junctions can reach an optimal rectification efficiency 
of ^ 0.8 28,29 and ~ 7.5 , 9,29 respectively. On the other hand, 
even though a simple NIN junction cannot provide heat rec¬ 
tification, it is possible to obtain ft opt as high as 3000 by re¬ 
alizing an asymmetric NININ chain with the central island 
coupled to the phonon bath by means of a cold finger. The 
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latter provides an efficient channel through which the diode 
can release energy in the non-transmissive temperature bias 
configuration. 30 This design can be used also in a NINISIN 
chain, which offers the opportunity to boost the performance 
of a single NIS junction and to overcome the fabrication com¬ 
plexity imposed by the high resistance asymmetry required in 
a NININ chain. The superconducting energy gap plays the 
role of a thermal bottleneck at low temperatures, while for 
higher temperatures it can generate an inversion of the recti¬ 
fication regime, as experimentally proven in Ref. 31 . Finally, 
we analyzed a NINIS 2 IS 1 IN chain as a testing ground for the 
experimental demonstration of a phase-tunable thermal diode, 
providing a larger rectification efficiency and requiring lower 
bias temperatures with respect to those obtained in previ¬ 
ous theoretical proposals. 9,29 These thermal rectifiers could be 
easily implemented by standard nanofabrication techniques 31 
and, combined with heat current interferometers, 7,8 might be¬ 
come building blocks of coherent caloritronic nanocircuits. 7-9 
Moreover, they could be connected to electronic coolers and 
radiation detectors, 1 2 3 4 5 6 7 8 9 10 offering the possibility to evacuate dissi¬ 
pated power in an unidirectional way. Finally, this technology 
might also have a potential impact in general-purpose cryo¬ 
genic electronic microcircuitry, e.g., solid-state quantum in¬ 
formation architectures. 11 12 13 
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